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ABSTRACT 

We present a large sample of Ly-a emitting galaxies spectroscopically con- 
firmed at redshift z ~ 4.5, based on IMACS spectroscopic observations of can- 
didate z ~ 4.5 Lya-emitting galaxies in the Large Area Lyman Alpha (LALA) 
narrow band imaging survey Cetus field. We identify 110 of them as z ~ 4.5 Lyo; 
emitters based on single line detections with no continuum emission blueward of 
the line. Six foreground galaxies are identified, either based on multiple lines or 
blueward continuum emission. The Ly-a confirmation rate varies from <50% to 
76% for candidates selected in different narrowband filters at slightly different 
redshifts. We find a drop in the LAE density at redshift 4.50±0.03 from redshift 
4.39 - 4.47 by a factor of 66%, which could be a large scale void in the distri- 
bution of star-forming galaxies (~ 18 Mpc along the line of sight and ~ 80 Mpc 
across). The sample includes many objects with equivalent widths (EW) > 200A. 
These large EW candidates are spectroscopically confirmed at the same rate as 
candidates with more modest EWs. A composite spectrum of all 110 confirmed 
LAEs shows the characteristic asymmetry of the Ly-a line. It also places new 
and stringent upper limits on the CIV 1549A/Ly-a and Hell 1640A/Ly-a line 
ratios, providing a new upper limit on the fraction of active galactic nuclei in 
Ly-a selected galaxy samples, and on the contribution of Pop III populations. 
Finally, we calculate the Ly-a luminosity function for our z ~ 4.5 sample, which 
is consistent with those at other redshifts, showing that there is no evolution in 
Lya luminosity function from z — 3.1 - 6.6. 
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1. Introduction 

The number of known high redshift galaxies has grown rapidly in the past decade, thanks 
to efficient wide field cameras and effective selection techniques (i.e., the dropout technique 
and Lya line search technique). Narrowband imaging surveys have become the most efficient 
approach for finding high redshift Ly-a emission lines galaxies. Great successes have been 
achieved recently in detecting high redshift Lya emitters using this technique (e.g. from z = 
3 to z = 6.96, see lye et al. 2006; Rhoads et al. 2000, 2003, 2004; Wang, Malhotra & Rhoads 
2005; Tapken et al. 2006; Westra et al. 2006; Shimasaku et al. 2004, 2006; Kashikawa et al. 
2006; Ajiki et al. 2006; Ouchi et al. 2005; Hu et al. 2002, 2004; Dawson et al. 2004, 2007; 
and references therein). Large well defined samples of high redshift Lya emitting galaxies 
are essential to study the nature of these galaxies (stellar population, mass etc, e.g. Lai et al. 
2008; Pirzkal et al. 2007; Finkelstein et al. 2007; 2009a), the mechanism of Lya production 
(Malhotra & Rhoads 2002) and escape (Finkelstein et al. 2008), the luminosity function and 
its evolution (Dawson et al. 2007; Malhotra & Rhoads 2004, 2006; Ouchi et al. 2008), and 
spatial clustering of star forming galaxies in the early universe (Wang et al 2005; Ouchi et 
al 2005; Malhotra et al. 2005; Kovac et al. 2007). 

The Large Area Lyman Alpha (LALA) survey (e.g., Rhoads et al. 2000) has obtained 
deep narrowband images in several fields, yielding a large sample of Lya emitters at z ~ 4.5 
(Dawson et al. 2004; 2007), and z « 5.7 (Rhoads & Malhotra 2001; Rhoads et al. 2003; 
Wang et al. 2005), plus a smaller sample at z ~ 6.5 (Rhoads et al. 2004). In this paper, 
we present the spectra of 110 confirmed Lya emitters at z ~ 4.5 in LALA Cetus field with 
the Inamori-Magellan Areal Camera and Spectrograph (IMACS, Dressier et al. 2006) on 
the 6.5 meter Magellan I Baade Telescope. Together with our previously confirmed sources 
from Keck observations (Dawson et al. 2004, 2007), we now have an overall sample of 171 
spectroscopically confirmed Lya emitters at z ~ 4.5. 

2. Observations 

2.1. Candidate Selection through Narrow-Band Imaging 

With the 8192x8192 pixel Mosaic CCD camera on the 4m Mayall Telescope at Kitt Peak 
National Observatory, we have obtained five deep narrow band images (36' x 36') in the LALA 
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Cetus field (02:05:20 -04:55; J2000.0) to search for candidate Lya-emitters at z « 4.5. The 
five narrowband filters utilized are overlapping with full width at half maximum (FWHM) 
80A (see Fig. 1 of Dawson et al. 2004). The central wavelengths are 6559, 6611, 6650, 6692, 
and 6730A for narrowband filter HaO, Ha4, Ha8, Hal2 and Hal6 respectively. The filters 
thus cover a total redshift range of 4.36 < z < 4.57 for Ly-cn, with approximately 50% overlap 
in bandpass between adjacent filters. The average seeing in the final stacked narrowband 
images are 1.00", 1.26", 1.12", 1.27" and 0.99", and the 5 a (2.4"diameter aperture) limiting 
magnitudes are 24.7, 24.8, 24.8, 24.7 and 24.1 (Vega) respectively. The zero points of the 
narrowband images were obtained by matching the narrowband photometry to that of the 
underlying broadband (R) image. 

To study the continuum emission properties of our candidates, we used deep broadband 
images in B w , R, I band from the NOAO Deep Wide Field Survey (NDWFS; Jannuzi 
& Dey 1999). We followed reasonably standard practices in our narrowband imaging data 
reduction and our Ly-a galaxy candidate selection, which are described in detail elsewhere 
(Rhoads et al 2000, Rhoads & Malhotra 2001, Wang, Rhoads & Malhotra 2005). Briefly, 
our selection criteria are: A) > 5a detection in a narrowband filter; B) a narrowband excess 
of > 0.75 magnitude, so that > 50% of the narrowband flux comes from an emission line; 
C) significance of the narrowband excess > 4a; and D) < 2a detection in By/ band. 

We identified a total of 226 candidate Lya-emitters at z ~ 4.5 in our primary sam- 
ple. The numbers of candidates selected in each narrow band filter are presented in Table 
1. Additionally, we performed an extended catalog selection using the spectral overlap of 
adjacent filters to go to fainter flux levels. In this case, we extended the first criterion to 
> 3(7 detection in one narrowband filter but > 5a when stacking two spectrally overlapping 
filters. Other criteria are the same. This leads to another 38 candidates, which comprise our 
extended catalog. 

2.2. Spectroscopic Observations 

In 2005 September 3-5 we obtained spectroscopy of 194 candidate z ~ 4.5 Ly-a emitters 
in the LALA Cetus field (167 from our primary sample, plus 27 from the extended catalog) 
using the Inamori-Magellan Areal Camera and Spectrograph (IMACS, Dressier et al. 2006) 
on the 6.5 meter Magellan I Baade Telescope. We used the short camera (f/2, with a 27.2' 
diameter field of view), and used the 300 line grism (Xuaze = 6700 A, 1.34 A/pixel) in 
combination with a 5650-9200A blocking filter to prevent overlap of higher order spectra. 
All observations employed 10.0" x 1.0" slitlets. The IMACS f/2 pixel scale is 0.2"/pixel. The 
corresponding spectral resolution is ~ 7A (though it could be somewhat better for compacat 
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objects observed in subarcsecond seeing). Five multislit masks (see Table 2) were observed 
for 2.5 to 3.3 hours in 0.5 hr increments. Thanks to the large field of view of IMACS, we 
were able to target around 50 candidate Lya-emitters in each mask, along with roughly 120 
other targets. Of these, 45 candidates were covered by more than one masks. 

The data were reduced using the IMACS version of the Carnegie Observatories System 
for MultiObject Spectroscopy (COSMOS) data reduction packagqj. We first determined 2d 
wavelength solutions for each science exposure using arc lamp exposures taken immediately 
before or after. We found relatively large wavelength residuals (~ 6 pixels) in the course of 
this calibration, perhaps pointing to inaccuracies in the geometric coefficients describing the 
layout of the detector chips in the focal plane. Regardless of their origin, these residuals are 
repeatable from exposure to exposure. They therefore do not significantly affect our search 
for Lya emission lines, though they do lead to extra uncertainties in the line wavelengths we 
determine (see below for details) and sky substraction. 

After wavelength calibration, each frame was first bias-subtracted and flat fielded. We 
then performed sky subtraction following the algorithm described by Kelson (2003), and 
extracted ID spectra from the 2D spectra using the task "extract-2dspec" for each slit. 
Frequently, this procedure produced an apparently "negative" blue continuum in the ID 
spectra of our confirmed Lya emitters, indicating small systematic errors in the software's 
measured sky level. In these cases, we checked that the continuum level appears to be zero 
in the 2D spectra, and then applied an additive correction to the ID spectrum. We note 
that none of our subsequent scientific results rely on the continuum level measured in the 
spectroscopic data. 

To control for possible spatial shifts along the slits between individual exposures of a 
mask, we measured the trace locations of the 30 brightest continuum sources in the mask. 
We corrected for any measured shifts while stacking the exposures for each mask, to increase 
the quality of the stacked 2D spectra. We also identified and removed cosmic ray hits while 
stacking the multiple exposures of each mask. Finally, we extracted the ID spectra (2A per 
bin) using the gaussian extraction algorithm. 



1 http://llama.lco.cl/~aoemler/COSMOS.html 
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3. Spectroscopic Results 

3.1. Line Identification 

We identified 110 single emission lines as Lya at z ~ 4.5, among our 194 observed 
candidates. Of these, 97 are new spectroscopic confirmations. We present a catalog of the 
110 confirmed LAEs at z ~ 4.5 in LALA Cetus field in Table 3. In Fig. [I] we presented the 
narrowband flux distributions of the targeted and confirmed LAEs in different narrowbands, 
and in Fig. [2] we present a set of sample IMACS spectra of confirmed LAEs. 

We also identify 6 foreground sources which either show multiple lines or blue continuum. 
The remaining 78 targets were classified as nondetections. Of the nondetections, there are 10 
whose expected emission line wavelength range (based on the narrowband data) is covered 
or partially covered by the IMACS ccd chip gap (which has a width of ~ 80 A). These 
ten include one source previously confirmed as a Lya emitter with Keck LRIS spectrum by 
Dawson et al. (2004). The non-detections could 

Following Rhoads et al. (2003, 2004) we perform two measurements of the line asym- 
metry. For both measurements, we first determine the wavelength of the emission line peak 
(A p ), and the data points where the line flux exceeds 10% of the peak on the blue side (Aio,&) 
and the red side (Aio, r )- The "wavelength ratio" is defined as a\ = (Aio,r-Ap)/(A p -Aio,b), and 
the "flux ratio" a/ as the ratio of the accumulated line flux between Aio, r and X p to that of 
between X p and Aio,6- In Fig. [3] we plot a\ versus aj for the 110 single lines emitters. The 
errorbars are obtained through one thousand Monte Carlo simulations in which we added 
random noise to each data bin. Both measurements are expected to be > 1 for asymmet- 
ric profile of high redshift Lya emission line due to the absorption by neutral hydrogen in 
outflows from star-forming galaxies and IGM. 

Due to the low resolution and limited spectral quality of our IMACS spectra, for most 
of the lines, we are unable to secure identification of asymmetric profile, i.e., with a a or a/ 
>1 at statistically significant level. This can be clearly seen from the large errorbars of a\ 
and cif in Fig. [31 Considering the large errorbars, we find all but 1 sources statistically 
satisfy (within la) ax > 1 or a,f > 1, and all but 3 sources satisfy both. This indicates that 
the spectral quality is insufficient for us to secure or exclude Lya identification based on line 
asymmetry measurements. 

Possible contaminations to our identifications of Lya lines are [O II] A3727 and [O III] 
A5007. While [O III] A5007 could be typically identified by neighboring [O III] A4959 and 
H/3, [O II] A3727 could be identified by looking for corresponding H/3 and [O III] A5007 at 
longer wavelength. These identifications have been done for more than 90% of the spectra. 
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For the rest the search for H/3 and [O III] A5007 is not possible because of the limited spectral 
coverage (or due to strong sky emission lines). For these sources, secure identification of [O 
II] A3727 required higher spectral resolution to resolve the doublet, which is not possible by 
our IMACS spectra. Luckily, most of the foreground [O II] A3727 emission line galaxies have 
been excluded from our candidate list by excluding sources detected in blueward broadband, 
and the contamination of [O II] A3727 in spectroscopically securely confirmed sample was 
verified to be very low (e.g. Dawson et al. 2004; 2007). For instance, with Keck DEIMOS 
spectra, Dawson et al. (2007) have identified two [O II] A3727 versus 59 Lya lines based on 
line profile, among candidate LAEs mostly selected similarly in LALA Bootes field. In this 
paper, we identify all the 110 single lines as Lya lines. 



3.2. Comparison with Dawson et al. 

With previous Keck observations, Dawson et al. (2004; 2007) have confirmed a total of 
17 LAEs at z ~ 4.5 in LALA Cetus field. While designing our IMACS masks, we gave lower 
priorities to these previous confirmed LAEs, in order to maximize the total number of LAE 
candidates with spectroscopic coverage. Thanks to the large field of view of IMACS, we still 
have 13 of previous confirmed LAEs covered by IMACS observations. 12 of the 13 sources 
are confirmed at z ~ 4.5 independently with the IMACS observations. The only one without 
IMACS confirmation has the expected line wavelength covered by IMACS chip gap. In Fig. 
H] we plot the observed Keck redshifts versus IMACS redshifts for these 12 sources, which 
demonstrates the uncertainties in the IMACS wavelength calibration {bz\ a = 0.005 at z pa 
4.5). Our IMACS exposures also covered 7 candidates in Cetus field which were previously 
reported as non-detections with Keck spectra (Dawson et al. 2004; 2007). None of these 
candidates showed any emission line in the IMACS spectra either. These results suggest that 
our IMACS spectra achieve sensitivities comparable to our previous Keck observations. The 
17 confirmations reported by Dawson et al (2004; 2007) were drawn from a total of 26 LAE 
candidates observed with Keck in the LALA Cetus field, for a success rate of 65%. We note 
that 24 out of these 26 candidates were selected from HaO, Ha4 or Ha8 band (since final, 
full-depth reductions of the Hal2 and Hal6 filter images were not available at the time of 
the Keck observations). For the HaO, Ha4 and Ha8 bands alone, the success rate of Dawson 
et al. is 71%. Our IMACS success rate for these same three filters is 66%, consistent with 
Dawson et al. 
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3.3. Equivalent width distribution 

In our previous studies among others(e.g. Malhotra & Rhoads 2002; Dawson et al. 
2004; 2007; Kudritzki et al. 2000), large fractions of candidate LAEs were reported to show 
rest-frame Lya line EW above 240 A, which is too large for normal stellar populations, 
unless they have a top-heavy initial mass function (IMF), zero (or very low) metallicity, 
and/or extreme youth (age < 10 6 yr). Significant contributions of type 2 active galactic 
nuclei (AGNs) to the large line EW has been ruled out (see §4.2 below). The large line 
EW could also likely be due to the clumpy and dusty inter-stellar medium in the galaxy 
which could attenuate the continuum photons but resonantly scatter the line photons at the 
surface of the dusty clouds, thus increase the observed EW (Neufeld 1991; Hansen & Oh 
2006; Finkelstein et al. 2007; 2008; 2009). 

Could the candidate LAEs with higher EW based on photometry (i.e. these with no 
continuum detection) be contaminations due to noise spikes, or spurious features? Our large 
spectroscopic sample enable us to perform a statistical test on this issue. In Fig. \5\ we 
plot the distributions of the broad to narrow band flux density ratio for our confirmed z ~ 
4.5 LAEs, foreground sources and spectroscopically undetected sources. The Lya line EW 
is a monotonic decreasing function of the broad to narrow band flux density ratio. In the 
figure we can see that these spectroscopically confirmed LAEs and undetected sources have a 
consistent distributions of broad to narrow band flux density ratio. A Kolmogorov-Smirnov 
test shows that two distributions are indistinguishable (p = 85%). This demonstrates that 
high EW candidates are confirmed at the same rate. 

In table 3 we also present the line equivalent width for all the spectroscopically con- 
firmed z fa 4.5 galaxies. Due to the large uncertainty in estimating the continuum level 
from optical spectra, we utilize narrow and broad band photometry to derive the line equiv- 
alent width. Following Maholtra & Rhoads (2002), we calculate the rest frame line EW 
as F Lya /f x /(l+z) = (W R N - W N R)/(R - N)/(l+z), where F Lya is the Lya line flux, f x the 
continuum flux density, Wjv and Wr the widths of the R band (1568 A) and the narrow band 
(80A), and R and N the integrated fluxes in the R band and narrow band respectively. 

Consistent with Malhotra & Rhoads (2002), we find only 50% of the confirmed LAEs 
show measured line EW < 240A in the rest frame. The rest either show larger EW or no 
detection of continuum at all. It's clear that the uncertainty in line EW mainly comes from 
that of the continuum, especially for high EW sources, we thus calculate the uncertainty of 
the line EW by applying la errorbar to the measured R band flux. In many cases that the 
upper limits of the line EW can not be constrained, lower limits are given. 

We point out that the line EW distributions of different selected samples are sensitive 
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to the selection criteria adopted. Particularly, the requirement of broad band detection (i.e. 
Hu et al. 2004) would exclude very high EW sources. For instance, only 1 out of 7 LAEs at 
redshift 3.7 selected by Fujita et al. (2003) in Subaru/XMM- Newton Deep Field show rest 
frame line EW above 200A (based on narrow and broad band photometry), however, this is 
after they excluded the candidates without broadband detection, most of which (6 out 7), 
if located at redshift 3.7, show line EW above 200A, consistent with LALA surveys. The 
EW distribution is also likely sensitive to the relative depths of the broad and narrow band 
images, e.g., deeper broad band images could enable detection of more sources with smaller 
EWs because more sources with stronger continuum could pass the selection due to smaller 
broadband flux uncertainty. 

We have presented above that most of the sources with high EWs have no broadband 
detection at all. It's interesting to note that if we sum the broad and narrow band fluxes for 
all the confirmed LAEs, we obtained an average EW of 182^1. If we sum the broad and 
narrow band fluxes for the confirmed LAEs with < 3cr continuum detection (100 out of 110), 
we derived an average EW of 316^g2 2 A. This also demonstrates that a large population of 
these individual LAEs have high EW. 

One final concern on the weak continuum level is whether they are physically weak or due 
to selection bias. That is, did the selection favor the sources whose broad band fluxes were 
underestimated due to random fluctuations? If this were an important effect, there should 
be one immediate consequence: these high EW sources would be systematically brighter in 
an independent broadband image. We adopted the independent broadband images obtained 
with MMT (Finkelstein et al. 2007; Wang et al. 2007) to test this possibility. In Fig. M we 
plot CTIO R band flux versus MMT r' "flux" for 44 confirmed LAEs with r' band coverage. 
Here the MMT r' "flux" was scaled by a factor which was derived by requiring that the 
clipped average R — r' color be zero for the full list of sources well detected in both images. 
We note that the band pass of MMT r' band is not identical to that of CTIO band, however, 
considering MMT r' band has much better efficiency between 6500A to 7200A than CTIO R 
band, while CTIO R band extends to 8000A with transmission below 40%, the fluxes from 
two bands are comparable for our purpose. From the figure we can see that these LAEs 
have consistent r' band fluxes, suggesting the possible selection bias we mentioned above, if 
there is any, is weak. The clearest exception is that fluxes that appear negative in the R 
band image (and whose R fluxes are therefore surely underestimated due to photon noise) 
mostly have r' fluxes that are zero or slightly positive. EW calculations show that 24 out 
of 44 confirmed LAEs with r' band flux measurements have EW above 240A based on old 
R band measurements, while 21 of the 44 show EW above 240A based on r' band fluxes. 
The EW calculations based on r' band fluxes have been adjusted to the r' band transmission 
curve. We conclude that the high EWs we have obtained are physically real, and cannot be 
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fully ascribed to selection bias. 

3.4. The stacked spetrum 

Stacking a large number of spectra could further help us understand the nature of the 
line emitters, especially for those with low signal to noise ratio. In Fig. [7J we present the 
stacked spectrum of all 110 confirmed Lya emitters at z ~ 4.5. We perform a variance- 
weighted co-addition of the 110 spectra. In each wavelength bin, data points from each 
spectrum (normalized by the peak flux of the Lya 1216A emission line) were averaged using 
two sigma clipping algorithm (one iteration) to remove artifical features due to sky line 
residuals, ccd edges, etc. 

In the composite spectrum, the only visible line feature is the asymmetric Lya line 
(highlighted in the cutout of Fig. [7]). In the co- added spectrum, we obtained a\ = 1.46±0.16, 
and (if = 1.42±0.09, consistent with those expected from the asymmetric profile of high 
redshift Lya emission line. The errorbars on a\ and a/ are also estimated based on large 
number of Monte Carlo simulations in which we added random noise to each data bin. 

In the co-added spectrum, we searched for possible CIV A1549 and He II A1640 line in 
the vicinity of expected wavelengths, however, neither of them could be detected at above 
2a level. Based on the flux uncertainties in the spectra, we obtained 2a detection limit (by 
fixing the line width to that of Lya) to the flux ratio of /(CIV A1549)//(Lya) and /(He II 
A1640)//(Lya), which are 3.7% and 7.4% respectively. These constraints are significantly 
stronger than we previously obtained based on 11 LAEs spectra with Keck (8% and 13% 
respectively, Dawson et al. 2004). 

4. Discussion 

4.1. Large scale structures 

We find significantly different success rates and total numbers of candidates obtained 
in images from different narrowband filters (see Table 1). The small number of candidates 
in the Hal6 filter can be ascribed to the shallower depth in the Hal6 narrow band image, 
and the lower success rate could thus be due to a higher percentage of contamination from 
noise spikes (assuming the total number of spurious detections due to noise spikes per image 
remains constant) and asteroids (with a small number of frames, asteroids could leave more 
spurious features in the stacked image). 
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The Hal2 narrowband image has the same depth as HaO, Ha4 and Ha8 bands, yet 
has many fewer candidates detected. This indicates that the source density at the redshift 
4.50±0.03 is much smaller than that in the redshift range of 4.39 - 4.47 (covered by HaO, Ha4 
and Ha8 filter). The success rates and number of candidates in HaO, Ha4 and Ha8 filters are 
similar, with Ha8 having marginally more candidates, and a marginally lower confirmation 
rate. 

In Fig. [8] we plot the spatial distributions of the candidates selected in each narrowband 
image. Similar large scale structure has been reported at various redshifts elsewhere (e.g. 
see Ouchi et al. 2005, Wang et al. 2005, Palunas et al. 2004, Kovac et al. 2007, and 
more). Quantitative analysis on the large scale structures presented in Fig. [8] requires 
careful examination on the candidate selection effect, the coverage of mask design, and the 
success rate of each mask in each narrow band, etc, which is beyond the scope of this paper. 
However we note that the amplitude of the drop in LAE source density at redshift 4.50±0.03 



(66%)g is consistent with the cosmic variance calculation (Trenti & Stiavelli 2008), which 
gives a 40% cosmic variance (la) for the volume of our survey. 



Significant contributions of AGNs to the large line EW has been ruled out with deep 
X-ray exposure and optical spectra (e.g. see Malhotra et al. 2003; Wang et al. 2004; Dawson 
et al. 2004; Gawiser et al. 2006). Deep X-ray images can be used to test whether these LAEs 
are AGNs, since they are expected to be strong X-ray emitters based on their Lja emission 
line. However, with two deep Chandra ACIS exposures (~ 180 ks each), we were unable 
to detect any of the 101 LAE candidates in LALA Bootes and Cetus field. The stacking 
technique helps to establish a strong constraint (3a) to the average X-ray luminosity of 
< 2.8 x 10 42 ergs s _1 , and we estimated that < 4.8% of the LAEs in LALA Bootes 
and Cetus field could be possible AGNs (Wang et al. 2004). No AGN lines have been 
identified in the obtained individual and composite spectra either (Dawson et al. 2004; and 
this paper). Dawson et al. (2004) have presented an upper limit of CIV to Lya line flux 
ratio of 8% (2a) based on the composite spectrum of 11 LAEs at z ~ 4.5. 

In §3.4 we have presented the composite spectra of 110 Lya emitters at z ~ 4.5. Even in 



2 We first calculated the expected numbers of LAEs in HaO, Ha4, Ha8 and Hal2 band respectively, by 
multiplying the total number of candidates with the spectroscopical success rate in each band. The fraction 
66% was obtained by comparing the expected number of LAEs in Hal 2 band with the average value in HaO, 
Ha4 and Ha8 band. 




4.2. Active Galactic Nucleus(AGN) Fraction 
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the composite spectrum, CIV line is invisible. A much stronger upper limit to the CIV line 
flux is derived, which is < 3.7% (2a) of the Lya flux, while typical ratio of CIV to Lya line 
in type II AGN^I is 22% (Ferland & Osterbrock 1986). This confirms previous results that 
the AGN contribution to LAEs is rather small. The formal upper limit from the CIV line 
is < 3.7/22 = 17%. This is somewhat weaker than our earlier X-ray based limits (Malhotra 
et al 2003, Wang et al 2004), but is an entirely independent constraint. Interestingly, this 
limit is below the AGN detection fraction for nearby (redshift z ~ 0.3) LAEs, reinforcing 
the observation of AGN fraction evolution by Finkelstein et al (2009b). 

4.3. Population III stars 

The search for the first generation (or Pop III stars) is one of the major targets in 
observational studies. However, until now there is no strong evidence of the detection of Pop 
III stars, and the search for such objects remains a major observational challenge (Schaerer 
2008). Recent numerical calculations have shown that since Lyman continuum flux increases 
with decreasing metallicity in stellar populations, very low metallicity stellar population can 
produce strong Lya line with maximum EW > 400 - 850A for Z < 10~ 5 ZQ (e.g. Schaerer 
2003; Tumlinson, Shull & Venkatesan 2003). Besides strong Lya emission line, Pop III or 
very low metallicity stellar population is also expected to produce strong He II A1640 due 
to its much stronger He II ionizing continuum. Thus many studies have been performed 
to search for strong He II A 1640 emission line in the spectra of high redshift LAEs, but no 
positive detection has been obtained yet. 

Consistent with Malhotra & Rhoads (2002) and Ouchi et al. (2008), we have found a 
large population of spectroscopically confirmed LAEs with rest frame Lya line EW above 
240A. We have stated above that these high EWs are real. This is in contrast to low EWs 
found at z = 3.1 & 2.2 (Gronwall et al. 2007; Nilsson et al. 2009), whether this is due to 
cosmic evolution or different selection criteria is not clear. We have also ruled out significant 
contribution of AGN to the large EW (§4.2). Another possibility to the observed large Lya 
EW is due to very metal-poor or Pop III stars. However, in both the individual and the co- 
added IMACS spectra, we detect no evidence of He II A1640 emission line. In the co-added 
spectrum, the He II A1640 emission line flux is constrained to be < 7.4% (2a) of that of 
Lya. Such an upper limit can be converted to log (Qne+fQii) < - 1.4, where Q,He+ and Qh 
are ionizing fluxes for He+ and H respectively (Schaerer 2003). This upper limit is close to 



3 Type 1 AGN could be ruled out in most (if not all) LAEs because of the lack of broad emission lines in 
the spectra. 
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the maximum predicted values for metal-free populations with IMFs including very massive 
stars (up to 500 M©, see Fig. 5 of Schaerer 2003). This indicates that although He II A1640 
is not detected in our composite spectrum, the derived upper limit is insufficient yet to rule 
out the existence of Pop III (or very metal-poor) stars in our LAEs at z ~ 4.5. Considering 
the fraction of Pop III galaxies among LAE sample at z ~ 4.5 could be rather small (such 
as a few percent, e.g. Scannapieco, Schneider & Ferrara 2003), the derived upper limit can 
not even rule out some model which predict comparable He II A1640 and Lya line fluxes for 
Pop III populations (e.g. Tumlinson et al. 2003). 



We have presented a large spectroscopically confirmed LAE sample at z ~ 4.5. In 
this section we present the Lya luminosity function based on this sample. Due to the 
large uncertainty in flux calibration in spectroscopic observations, we adopt photometric 
measurements to calculate the Lya fluxes. The Lya fluxes were derived from narrowband 
and broadband photometry, where f Lya = f nb - f continuum* w nb and f continuum = {firfnb)/ (wr- 
w n b). The line luminosities of the confirmed LAEs were calculated assuming a cosmology 
model of HO = 71 km/s/Mpc, fl m = 0.27 and fi A = 0.73. 

For each confirmed LAE, we first calculated the comoving volume V max over which the 
source could be selected by our survey from z min to z max , which are 4.36 and 4.65 respec- 
tively through band HaO to band Hal6. The survey area is 35' x 35'. Effects of different 
flux limits in narrow bands were also taken into account. We show the incompleteness of 
our spectroscopical observations in Fig. In the figure we plot the fraction of spectroscop- 
ically targeted LAEs as a function of Lya line flux. Such incompleteness was taken into 
account while calculating the luminosity function. To be conservative, we assume that all 
the spectroscopically non-detected sources are not LAEs. 

The derived Lyct luminosity function of LAEs at z ~ 4.5 is plotted in Fig. [TU1 Following 
Malhotra & Rhoads (2004), we fit the luminosity function with a Schechter function 



During the fit, we fixed the slope a at -1.5 since our data is insufficient to have it constrained. 
The best-fit parameters are L* = 6.3 ± 1.5 x 10 42 erg s" 1 and $* = 3.5 ± 1.3 x 10~ 4 Mpc~ 3 . 
Considering the large uncertainties in the measurements and differences in Lyct emitting 
galaxies selection criteria among different research groups, the best-fitted L* and $* are 
consistent with the range of values found in various surveys at 3 < z < 6.5 (see the lower 
panel in Fig. [TO]) . 



4.4. Lya Luminosity Function 




(1) 
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5. Summary 

We have presented the largest spectroscopic sample of z ~ 4.5 Ly-a emitting galaxies 
yet studied. Examining the properties of this sample, we reach several conclusions: 

(1) Large equivalent width Ly-a emitters, which we have earlier reported as a common 
feature of narrowband-selected samples (Malhotra & Rhoads 2002), are spectroscopically 
confirmed at the same rate as less spectacular line emitters. Thus, the high equivalent width 
sample is not dominated by noise spikes. Nor can the large equivalent widths be solely a 
byproduct of calculating the EW using the same broad band data employed in candidate 
selection, because our large EW subsample changes little when we calculate EW using an 
independent broad band image completely independent of object selection. 

(2) The CIV 1549A line is (on average) < 3.7% of the Ly-a line strength, placing 
an upper limit of < 17% on the AGN fraction in the Ly-a selected sample. This limit is 
independent of the X-ray to Ly-a ratio in AGN and therefore complements earlier limits on 
the AGN fraction based on X-ray photometry. 

(3) The Hell 1640A line is < 7.4% of the Ly-a line strength. This is the most stringent 
limit to date on the Hell emission from these objects. Though it does not strongly rule out 
many Population III models yet, it is nearly deep enough to do so and could be combined 
with a few similarly sized samples in future. 

(4) We present a new Ly-a luminosity function, based on a large spectroscopic data 
set at z ~ 4.5. Allowing for differences in Ly-a galaxy selection criteria among different 
research groups, this luminosity function remains broadly consistent with others over the 
range 3 < z < 6.5. 

The work of JXW is supported by National Basic Research Program of China (973 
program, Grant No. 2007CB815404), Knowledge Innovation Program of CAS (Grant No. 
KJCX2-YW-T05), and Chinese National Science Foundation (Grant No. 10825312). Our 
data were obtained using community access telescope time made available under the National 
Science Foundation's Telescope System Instrumentation Program (TSIP), awarded by the 
National Optical Astronomy Observatory. The work of JER and SM is supported in part by 
NSF grant AST-0808165. 
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Table 1. Numbers of candidate, targeted and confirmed Lya-emitters at z 4.5 in each 
narrow band filter in the LALA Cetus field. The "ext" line represents candidates in the 
extended sample, which goes to a fainter flux level by combining data from filters that 

overlap in wavelength. 



filter 


average seeing 


limiting mag" 


candidates 


targeted 


confirmed 


HaO 


1.00" 


24.7 


57 


48 


35 


Ha4 


1.26" 


24.8 


64 


45 


34 


Ha8 


1.12" 


24.8 


80 


56 


33 


Hal2 


1.27" 


24.7 


38 


22 


9 


Hal6 


0.99" 


24.1 


31 


23 


8 


cxt 






38 


27 


12 



a 5 a (2. 4" diameter aperture) limiting magnitudes (Vega) 
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Table 2. Numbers of targeted and confirmed Lya-emitters at z ~ 4.5 in each mask 



mask 


exp time (ks) 


targeted 


confirmed 


1 


11.3 


52 


35 


2 


8.9 


58 


33 


3 


10.8 


44 


21 


4 


9.9 


58 


26 


6 


9.0 


46 


17 
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Table 3. IMACS Spectroscopically confirmed Lya-emitters at z m 4.5 in the LALA Cetus 

field 



Target z a Lya flux 6 EW^ sf (A) 



J020623. 


1-044851 


4, 


431 


2 


584±0 


.357 


>105 


J020616. 


,5-045457 


4, 


378 


3 


.795±0. 


410 


—58 


J020616. 


.2-050315 


4, 


428 


2 


399±0 


435 


70+^! 

— 33 


J020615. 


9-050555 


4, 


.380 


2 


,389±0 


353 


>82 


J020614. 


4-050523 


4 


.380 


7 


,448±0 


376 


' ^ — 126 


J020613. 


.1-045101 


4 


.372 


2 


648±0. 


.401 


>83 


J020612. 


.4-045643 


4, 


.385 


3 


029±0 


361 


>118 


J020556. 


.7-045854 


4, 


.421 


2 


,850±0 


.359 


xu ' —52 


J020556. 


.5-050917 


4, 


.388 


2 


,817±0. 


.408 


30+f 


J020544. 


.7-045921 


4, 


.382 


2 


,117±0. 


.380 


>78 


J020542. 


.1-050520 


4, 


.382 


1 


,993±0. 


377 


>77 


J020536. 


3-044822 


4, 


.414 


5 


.660±0 


399 


>179 


J020533. 


.0-050042 


4, 


395 


2 


.844±0 


408 


>125 


J020526. 


.2-050724 


4, 


390 


2 


,576±0 


377 


>67 


J020525. 


.2-050048 


4, 


.400 


4 


,276±0. 


.431 


186+if 9 


J020524. 


.5-045241 


4, 


.390 


3 


.221±0. 


.420 


>135 


J020518. 


.1-045615 


4, 


.406 


4 


,074±0. 


.417 


>161 


J020513. 


.9-044638 


4, 


395 


3 


.268±0 


.431 


>87 


J020510. 


.5-044456 


4, 


.406 


3 


513±0. 


434 


>70 


J020509. 


7-044005 


4, 


410 


5 


,030±0. 


372 


— 50 


J020500. 


1-050731 


4, 


388 


7. 


,290±0. 


.413 


83+ " :i 
— 18 


J020457. 


.3-051157 


4, 


.398 


5 


,928±0 


.414 


>116 


J020451. 


.6-045737 


4, 


365 


2 


.221±0 


.386 


>85 


J020450. 


.6-050107 


4, 


.378 


2 


.256±0. 


.399 


>87 


J020450. 


.9-044323 


4, 


.401 


2 


124±0 


.414 


49+"' 7 

— zU 


J020446. 


.9-050116 


4, 


.354 


2 


,208±0. 


.387 


>92 


J020438. 


.6-044437 


4, 


.398 


3 


,233±0. 


.428 


43+?^ 


J020434. 


6-051016 


4, 


.378 


2 


.813±0 


.410 


99+i" 


J020434. 


.5-050516 


4, 


401 


6 


.029±0. 


.452 




J020432. 


3-045519 


4, 


360 


2 


.740±0. 


.413 


>97 


J020428. 


.5-045924 


4, 


391 


2 


,766±0. 


.393 


yD -45 


J020427. 


.4-050045 


4, 


.390 


2 


486±0. 


.411 


>100 


J020425. 


7-045810 


4, 


383 


2 


.788±0. 


.400 


>68 


J020414. 


.3-051201 


4, 


.401 


4 


,139±0. 


.453 


>119 


J020413. 


.8-044703 


4, 


.385 


2 


.928±0. 


.422 


>72 


J020628. 


.0-045307 


4, 


.462 


2 


011±0 


.376 


>69 


J020623. 


8-044221 


4, 


444 


3 


.859±0 


.366 


>169 


J020603. 


5-050019 


4, 


441 


2 


.297±0. 


.357 


>59 


J020601. 


.7-050258 


4, 


.433 


3 


.134±0. 


.373 


>126 


J020554. 


,3-044147 


4, 


.428 


2 


.761±0 


.346 


>77 


J020550. 


6-045005 


4, 


429 


2 


.005±0 


.350 


69+3f 


J020536. 


.6-045520 


4, 


418 


2 


,689±0 


.393 


99 +b 5 f 


J020531. 


4-044053 


4, 


459 


2 


,115±0 


.360 


co + 130 
°°-23 


J020527. 


,7-043944 


4, 


.436 


3 


,848±0. 


.374 




J020519. 


.2-044902 


4, 


.421 


2 


.531±0 


.380 


>95 


J020517. 


.5-044347 


4, 


436 


2 


.420±0 


.396 


>99 
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Table 3 — Continued 



Target z a Lya flux 6 EW£ est (A) 



J020513. 


8- 


-043833 


4. 


.419 


2 


723±0 


376 


>117 


J020507. 


.2- 


-050636 


4 


.465 


2 


,861±0. 


.362 


>119 


J020506. 


8- 


-043858 


4. 


.429 


3 


873±0. 


355 


o-7+169 
87 -31 


J020458. 


0- 


-050140 


4. 


.452 


2 


.155±0 


377 


oo+411 
82 -42 


J020457. 


.3- 


-043847 


4. 


.431 


2 


086±0. 


.366 


>47 


J020452. 


1- 


-043802 


4. 


.431 


3 


,346±0 


361 


>88 


J020450. 


.2- 


-045828 


4. 


.423 


2 


510±0. 


366 


>94 


J020448. 


.0- 


-044257 


4 


.426 


2 


.871±0. 


.374 


7Q-f 200 
' y -33 


J020445. 


.5- 


-050008 


4. 


416 


2 


,297±0 


.383 


57+25 


J020439. 


.0- 


-045233 


4. 


.421 


2 


,387±0. 


.398 


>79 


J020438. 


.9- 


-044401 


4. 


.414 


2 


.151±0 


376 


>54 


J020438. 


9- 


-051116 


4. 


.447 


9 


920±0. 


366 


>247 


J020434. 


.9- 


-045703 


4. 


.433 


2 


,422±0. 


363 


. „_(_91 
48 -18 


J020432. 


3- 


-050916 


4. 


.456 


2 


992±0 


373 


>62 


J020429. 


1- 


-044307 


4. 


.434 


3 


992±0. 


.389 


seta! 3 


J020424. 


,7 


-043912 


4. 


.433 


3 


.839±0 


374 


>110 


J020418. 


.1- 


-050748 


4. 


.449 


2 


.813±0. 


.354 


>105 


J020617. 


.0- 


-045438 


4. 


.475 


2 


.081±0. 


366 


>65 


J020614. 


.0- 


-050031 


4. 


.479 


3 


,085±0 


354 


>131 


J020612. 


.3- 


-045608 


4. 


.467 


2 


,529±0. 


368 


>103 


J020611 


.1- 


-044601 


4 


.488 


2 


293±0. 


.349 


>84 


J020603. 


.8- 


-044905 


4. 


467 


3 


462±0. 


348 


47t™ 


J020557. 


,7 


-051003 


4. 


.457 


1 


901±0. 


.352 


>48 


J020557. 


.5- 


-045107 


4 


461 


2 


229±0. 


.334 


>85 


J020547. 


0- 


-045744 


4. 


.456 


2 


,785±0 


357 


>106 


J020536. 


.4- 


-050751 


4. 


484 


4 


419±0 


.339 


„_-f-192 
97 -35 


J020535. 


.6- 


-044430 


4. 


.498 


3 


,408±0 


413 


>142 


J020534. 


6- 


-044322 


4. 


507 


6 


,244±0 


391 


49 -9 


J020531. 


.2- 


-044704 


4. 


.459 


2 


,292±0 


.347 


>95 


J020526. 


.5- 


-050310 


4. 


469 


2 


,377±0. 


.346 


>82 


J020525. 


.7 


-045927 


4. 


.487 


1 


.962±0. 


.350 


>53 


J020525. 


.1- 


-044552 


4. 


.456 


2 


.688±0. 


351 


>91 


J020525. 


.1- 


-044515 


4. 


.444 


1 


778±0 


353 


>51 


J020456. 


6- 


-050145 


4. 


444 


1 


,796±0 


351 


46+H3 
* D -20 


J020454. 


.8- 


-045725 


4. 


.461 


2 


.146±0 


.348 


>47 


J020451. 


.3- 


■045147 


4. 


475 


2 


.000±0 


.342 


>44 


J020448. 


.3- 


-044628 


4. 


.474 


1 


,924±0. 


379 


>75 


J020442. 


.1- 


-050652 


4. 


.482 


2 


103±0. 


360 


>65 


J020432. 


.5- 


-044000 


4. 


470 


2 


153±0. 


.362 




J020429. 


.7 


-050251 


4. 


.462 


2 


152±0 


.377 


>70 


J020425. 


A- 


-045610 


4. 


461 


3 


,000±0. 


.361 


39l«? 


J020414. 


.6- 


-045956 


4. 


477 


2 


,777±0 


.352 


>104 


J020618. 


.2- 


-045024 


4. 


.523 


3 


,945±0. 


.406 


cc+84 
00 -15 


J020608. 


0- 


-044623 


4. 


512 


2 


,374±0 


.414 


>52 


J020604. 


.1- 


-050752 


4. 


495 


2 


.007±0 


.389 




J020557. 


6- 


-043614 


4. 


493 


3 


,290±0 


.605 


>81 
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Table 3 — Continued 



Target 




Lya flux 6 


EW= est (A) 


J020617.0-044833 


4.541 


3.583±0.679 


99±™ 


J020614.0-043652 


4.561 


4.533±0.796 


98+4J 3 


J020603.5-045637 


4.536 


4.303±0.659 


27+f 


J020602.0-045633 


4.559 


4.398±0.661 


136l« 9 


J020544. 1-045049 


4.543 


4.404±0.660 


>200 


J020503.1-044212 


4.515 


3.108±0.657 




J020534.7-051205 


4.423 


2.673±0.408 


>84 


J020501. 7-045639 


4.395 


1.729±0.378 


22+f 


J020442.7-044200 


4.408 


1.627±0.408 


45+ 114 


J020441. 3-044515 


4.424 


1.678±0.422 


>58 


J020427. 1-044736 


4.391 


2.028±0.404 


>77 


J020414.4-044411 


4.414 


1.845±0.425 


>55 


J020413.5-044040 


4.431 


2.049±0.376 


94tif 8 


J020618.2-043851 


4.442 


1.544±0.352 


qo+69 
dz -14 


J020612.3-045347 


4.442 


1.544±0.378 


>52 


J020512.9-045025 


4.439 


1.890±0.383 


>66 


J020513.2-045331 


4.484 


2.203±0.420 


>63 


J020606.9-045021 


4.539 


1.364±0.383 


>39 



a The redshift was measured from the wavelength of the 
peak pixel in the line profile from 1-d spectra. Due to the 
uncertainties in the wavelength calibration (see text for de- 
tails), we expect the error in the measurement to be 8z ~ 
0.005. 

b The line flux was measured based on narrowband and 
broadband photometry, in unit of 10 -17 erg.cm _2 .s _1 . 

c The EW were obtained based on broadband and narrow- 
band photometry. See §3.3 for details. 
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Fig. 1. — The narrowband flux distributions of IMACS targeted candidate (solid) and con- 
firmed Lja emitters (dashed). Different panels plot the distributions of candidates selected 
in different narrow bands. 
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Fig. 2. — Sample IMACS spectra for Lya line from the set of 110 confirmed LAEs at z ~ 
4.5 in LALA Cetus field. 
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Fig. 3. — The plot of two measurements of the line asymmetry for 110 Lya emitters (see 
text for definitions of a x and a/). Due to the low resolution of our IMACS spectra, both 
measurements suffer large uncertainties. Considering the large errorbars, all but 1 source 
satisfies (within la) a\ > 1 or a/ > 1, and all but 3 sources satify both. 
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Fig. 4. — The IMACS redshift versus Keck redshift for 13 previously confirmed LAEs with 
Keck spectra, demonstrating the uncertainty of IMACS wavelength calibration. 
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Fig. 5. — The distribution of the broad to narrow band flux density ratio for the spec- 
troscopically confirmed z » 4.5 sources, foreground sources, and undetected sources. The 
typical error in the flux density ratio is plotted in the bottom left corner. The line rest 
frame line EW (marked on the top of the plot) is a monotonic decreasing function of the 
flux density ratio. From the figure we see no difference in the flux density ratio distribution 
of pectroscopically confirmed z « 4.5 sources and undetected sources. The KS test yields 
no difference in these two populations (p = 85%). This suggests the EW distribution of all 
candidates or all confirmed z » 4.5 sources can represent that of all true z » 4.5 sources in 
our sample. 
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Fig. 6. — The flux density in the MMT r' band versus CTIO R band for 44 confirmed LAEs. 
Spectroscopical non-detections are also overplotted. The solid line represents equal lines 
in both bands. The typical errors in the measurements are presented in the bottom right 
corner. 
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Fig. 7. — The co-added spectrum of the 110 IMACS confirmed Lya emitters at z ~ 4.5 in the 
LALA Cetus field. Dashed lines show the one sigma uncertainty in the composite spectrum. 
In the co-added spectrum, neither CIV A1549 nor He II A1640 line could be detected in the 
vicinity of expected wavelengths, and 2a detection limits were given in the text. 
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Fig. 8. — The sky distribution of the z ~ 4.5 candidates (confirmed LAEs are over plotted 
as solid squares) detected in different narrowband. The size of the box is 38' x 38'. The 
distribution of the confirmed objects projected on the sky is subject to complex constraints 
on slit placement in multislit spectroscopy, so angular clustering in this figure should not 
be interpreted as a pure signature of galaxy clustering. Note except for Hal6 band, which 
is 0.6-0.7 magnitude shallower, all other four narrow band images have comparable limiting 
magnitudes. This figure clearly show that the source density of LAEs detected remains 
constant from redshift 4.39 to 4.47, but shows a significant drop at redshift 4.50. 
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Fig. 9. — The distribution of the narrow band fluxes for all candidate LAEs (solid line), and 
those targeted for spectroscopy (thick shaded region). For our narrowband filters, a flux 
of 1 microJansky corresponds to a line flux of 5.5xlCT 17 erg/cm 2 /s for an object with no 
continuum. The fraction of candidates LAEs with spectroscopic follow-up is over-plotted 
(dashed line). 
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Fig. 10. — The differential luminosity function of LAEs at z ~ 4.5. The best-fit Schechter 
function with slope fixed at -1.5 is over plotted (solid line). The best fit LF published in 
Dawson et al. (2007) is also over-plotted for comparison (dashed line). In the lower panel, 
we compare L* and $* from various surveys at different redshifts. W09: this work; D07: 
Dawson et al. (2007); O08: Ouchi et al. (2008); S06: Shamasaku et al. (2006); K06U: 
Kashikawa et al. (2006) upper limit; K06L: Kashikawa et al. (2006) lower limit. Best-fitting 
parameters obtained by fixing a = -1.5 are taken for all surveys except for D06 (a = -1.6). 



